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We have obtained deep optical images with the Very Large Telescope at ESO 
of the first well-localized short- duration gamma-ray burst, GRB 050509B. From 
V and R imaging, initiated ~ 2 days after the GRB trigger and lasting up to 
three weeks, we detect no variable object inside the small Swift/XBT X-ray error 
circle down to 2a limits of V = 26.5 and R = 25.1. The X-ray error circle includes 
a giant elliptical galaxy at z = 0.225, which has been proposed as the likely host 
of this GRB. Our limits indicate that if the GRB originated at z — 0.225, any 
supernova- like event accompanying the GRB would have to be over 100 times 
fainter than normal Type la SNe or Type Ic hypernovae, 5 times fainter than 
the faintest known la or Ic SNe, and fainter than the faintest known Type II 
SNe. Moreover, we use the optical limits to constrain the energetics of the GRB 
outflow. Simple models indicate that, unless the intrinsic energy in the outflow 
from GRB 050509B was <C 10 51 erg, there was very little radioactive material 
with efficient decay timescales for generating a large luminosity. These limits 
strongly constrain progenitor models for this short GRB. 

Subject headings: gamma rays: bursts — supernovae: general 



1. Introduction 

While it is now well established that long-duration 7-ray bursts (GRBs) coincide with 
the explosions of massive stars leading to very energetic core-collapse supernovae (SNe) 
(Galama et al. 1998; MacFadyen & Woosley 1999; Bloom et al. 1999; Stanek et al. 2003; 
Hjorth et al. 2003), the origin of short/hard GRBs, characterized as having durations < 2 s 
and hard spectra (Kouveliotou et al. 1993), remains unknown. There have been as yet no 
afterglow detections in the very few cases where searches for optical counterparts of short 
GRBs were performed, primarily due to the lack of early and precise localizations (Kehoe et 
al. 2001; Gorosabel et al. 2002; Hurley et al. 2002; Klotz, Boer & Atteia 2003). 

Recently, the Swift satellite (Gehrels et al. 2004) provided the first rapid and accurate 
X-ray localization of a short/hard GRB, opening the window for rapid progress on the 
origin of short GRBs. GRB 050509B (Gehrels et al. 2005) was detected on 2005 May 09 
at 04:00:19.23 (UT) by the Swift Burst Alert Telescope (BAT). It was a short (40 ms) and 



1 Based on observations collected at the European Southern Observatory, Paranal, Chile (ESO Programme 
075.D-0261). 
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fairly hard burst. The Swift X-ray Telescope (XRT) determined the source location to be 
R.A. = 12 h 36 m 13*58, decl. = +28° 59' 01"3 (J2000, error radius 973, Gehrels et al. 2005). 

The error region of GRB 050509B was observed by several groups (see, e.g., Bloom et 
al. 2005; Gehrels et al. 2005; Cenko et al. 2005; Castro-Tirado et al. 2005). Remarkably, 
the burst error circle overlaps with a giant elliptical galaxy, 2MASX J 1236 1286 + 2858580 
(hereafter Gl) at z = 0.225 (see Fig. 1 and Bloom et al. 2005), belonging to the cluster of 
galaxies NSC J123610+285901 (Gal et al. 2003). The probability of a chance alignment of a 
GRB and such a galaxy is of order ~ 10~ 3 . Assuming that Gl is therefore the host galaxy, 
Bloom et al. (2005) and Gehrels et al. (2005) argued that a likely origin of GRB 050509B is 
a neutron star (NS)-NS or NS-black hole (BH) merger. 

It should be noted that a merger does not necessarily imply the absence of optical or 
other long- wavelength phenomena after the GRB. For example, the 'mini-SN' model (Li & 
Paczyhski 1998) predicts a bright optical flash of much shorter duration than the one from 
a 'normal' SN, typically of about one day. But there are alternative scenarios for the origin 
of short GRBs. Zhang, Woosley & MacFadyen (2003) have suggested that short GRBs may 
be a variant of long GRBs, e.g., 'collapsar'-like events leading to stripped-core, core-collapse 
SNe, much like those seen in conjunction with long GRB afterglows (see also Ghirlanda, 
Ghisellini & Celotti 2004; Yamazaki, Ioka & Nakamura 2004). An alternative suggestion 
is that short GRBs are related to thermonuclear explosions, leading to Type la SNe (Dar 
& De Riijula 2004; Dado, Dar & De Rujula 2005). Finally, Germany et al. (2000) even 
suggested that the peculiar Type II SN 1997cy was related to the short GRB 970514 based 
on their temporal and spatial coincidence. It is obvious from the above, that a search for a 
SN associated with GRB 050509B would help constrain both the energetics of short GRBs 
and, possibly, their progenitor models (Fan et al. 2005). 

We have, therefore, obtained deep images of the XRT error circle at the expected peak 
time of the putative SN, as well as early images for comparison. In this Letter we present 
our observations and analysis (§ 2) and discuss the constraints these set on short GRB 
energetics and progenitor models (§3). A cosmology with H = TOkms^ 1 Mpc -1 , f2 m = 0.3, 
and £l A = 0.7 is assumed throughout this Letter. 

2. Observations and data analysis 

We obtained deep V and R band images containing the XRT error circle with the 
FORS instruments on the Kueyen and Antu 8.2-m unit telescopes of the ESO Very Large 
Telescope at several epochs (Table 1). The data obtained 8 and 12 days after burst were 
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strongly affected by the proximity of the Moon. Consequently our deepest images were 
obtained during the first and last sets of observations (a few days and ~3 weeks after burst). 

The data were bias subtracted, overscan corrected and flatfielded using morning skyflats. 
The photometric calibration was based on known FORS2 zero points and verified against the 
SDSS and photometry of the field obtained at the Tautenburg observatory. We estimated 
the field limiting magnitudes by doing photometry on a large (~ 50) number of objects. We 
used the IRAF task phot with an aperture of twice the seeing disk, and obtained the 3a 
limiting magnitudes given in Table 1 from the errors on the derived magnitudes. 

The first image, obtained 1.85 days after the burst, revealed a large number of very 
faint objects inside the XRT error circle, as well as Gl (see Fig. 1 and Hjorth et al. 2005; 
Gehrels et al. 2005; Bloom et al. 2005). The surface brightness of Gl within the error circle 
is however only at a level of 20% of the nightsky level in both the V and the R bands, with 
the exception of the central 1" of Gl, corresponding to less than 1 % of the error circle. The 
pixel-to-pixel photon noise is therefore only marginally affected by Gl. 

To search for variable objects we subtracted the early images from the late images in 
each band. The images were aligned, the sky background subtracted, and the images scaled 
to the same brightness level. We also convolved the image with the best seeing with a 
spatially variable kernel to match the inferior seeing of the other image, according to the 
method outlined in Alard & Lupton (1998). This provided very clean subtractions, except 
for near the center of the galaxy (Fig. 1). No variable sources were detected. 

The elliptical galaxy Gl has a simple morphology. It is therefore possible to fit and 
subtract a smooth model. This allows detection of superimposed point sources with a pixel- 
to-pixel noise which is reduced by a factor \/2, formally corresponding to a 0.38 mag deeper 
limit. A smooth fit of Gl was established by dividing the area around the galaxy into 
annuli of increasing widths, centered on the nucleus of the galaxy. A robust fitting technique 
was used to fit a harmonic series to pixel values within each of the annuli (Sodemann & 
Thomsen 1994). The full model was obtained by cubic spline interpolation between the 
harmonic coefficients in the radial direction. Finally the smooth model was subtracted from 
the galaxy image (Fig. 1). In the subtracted image we only detect one new faint (V ~ 26) 
object inside the XRT error circle, 0'.'78 east and 2'.'56 north of the Gl galaxy center. This 
object does not appear to be variable. No other object brighter than the already known 
sources (Bloom et al. 2005) is present inside the XRT error circle. 

To determine the detection limit in the images, a point-spread function was constructed 
from stars in the field. By superimposing artificial stars on Gl, we find that stars of V ~ 26.5 
and R ~ 25.1 are recovered in the difference images with 95% confidence, except within the 
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central 1" of the galaxy core. In the very nucleus of Gl, the galaxy subtraction is poor and 
we can only detect a source of V ~ 24. The corresponding limits for the galaxy subtracted 
images are 0.3-0.4 mag deeper. 

3. Discussion 

We plot the limits derived in § 2 in Fig. 2 along with a number of SN lightcurves as 
they would appear at z — 0.225. The Type la templates are from Nugent's compilation 2 and 
include (i) a template of a normal Type la SN (Nugent, Kim & Perlmutter 2002) and (ii) 
a template based on the very sub-luminous Type la supernovae SN 1991bg and SN 1999by. 
The Type Ic SNe plotted are (iii) the very energetic Type Ic SN 1998bw associated with 
the long GRB 980425 (Galama et al. 1998) and (iv) the faint, fast-rise Type Ic SN 19941 
(Richmond et al. 1996), which was not associated with a GRB but provides a good fit to 
the lightcurve bump in XRF 030723 (Fynbo et al. 2004). Figure 2 clearly demonstrates that 
even the faintest of these SNe would have been detected at the time of our observation at a 
level ~ 1.8 mag brighter than our limit (> 5.2 mag fainter than a SN like SN 1998bw). 

Type II SNe come in various flavors, the faintest of which are Type IIP. Our limit of 
V = 26.5 translates into M B > —13.3 at z — 0.225. All the SN peak magnitudes included 
in Richardson et al. (2002) are brighter than this, including the faintest Type IIP SNe. 

From the above we conclude that if GRB 050509B were associated with a normal bright 
SN, its host galaxy must either be at a high redshift (z > 1.2), consistent with the constraints 
on the redshifts of the faint galaxies in the XRT error circle (see Bloom et al. 2005) or, if 
it indeed is at z = 0.225, its SN light must have been extinguished by dust along the line 
of sight. The latter option appears unlikely as Gl is an elliptical galaxy with very little 
star formation (Bloom et al. 2005; Gehrels et al. 2005) and the likely background sources do 
not appear strongly reddened. We can therefore conclude that there was no SN of known 
type and characteristics associated with GRB 050509B if it occurred in Gl. There remains 
a (small) probability that GRB 050509B was at a much higher redshift. However, in the 
following we proceed under the assumption that Gl is the host galaxy. 

The absence of a SN rules out models predicting a normal SN la associated with short 
GRBs (Dar & De Riijula 2004; Dado, Dar & De Riijula 2005). Likewise, our observations dis- 
favour a GRB 050509B progenitor similar to long GRBs, i.e., a collapsar origin. Observations 
of long GRBs at z < 0.7 are consistent with all having SN bumps (Zeh, Klose & Hartmann 



2 http: / / supernova.lbl.gov/~nugent/nugent_templates.html 
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2004) and all GRBs at z < 0.4 have had SN features [GRB 980425, z = 0.0085, Galama et 
al. (1998); GRB 031203, z = 0.1055, Malesani et al. (2004); GRB 030329, z = 0.1685, Hjorth 
et al. (2003); GRB 011121, z = 0.362, Garnavich et al. (2003)]. The situation is less clear re- 
garding X-ray flashes (XRFs); a bright SN was associated with XRF 020903 (z = 0.251), but 
no SN (and no optical afterglow) was detected in XRF 040701 (with a probable z = 0.215) 
down to a limit at least 3 mag fainter than SN1998bw (Soderberg et al. 2005). 

We now proceed to use our derived limits to constrain the energetic properties of the 
outflow from GRB 050509B. Bloom et al. (2005) find that both the isotropic equivalent 
energy output in 7-rays, £ , 7) i SO , and the afterglow X-ray luminosity, Lx, of GRB 050509B are 
significantly smaller than those of long GRBs. The most straightforward conclusion is that, 
while GRB 050509B was a highly relativistic event (Gehrels et al. 2005), it was intrinsically 
less energetic compared to long GRBs, with relatively little energy [< 10 49 (f2/47r) erg] in 
highly relativistic ejecta with an initial Lorentz factor r > 100 (from E^) and with not 
much more energy in material with r > 3(£ , 5i/no) 1 ^ 8 , where £k,iso = 10 51 -E5i erg is the 
isotropic equivalent energy in the afterglow shock and n = tiq cm~ 3 is the external density. 3 
Moreover, the total observed energy from the burst was much smaller than the available 
energy in a NS-NS or NS-BH merger, or in most other progenitor models, suggesting that 
more energy was released in slower ejecta. The amount of energy in material above a certain 
initial four- velocity, E(> T p ), is very uncertain theoretically, but may be constrained by 
our late time upper limit on the optical emission. 

There are two ways to produce the most readily detectable emission from the outflow 
associated with GRB 050509B. It can either originate in the shock created by the outflow 
as it drives into the ambient medium, similar to both a long GRB afterglow for relativistic 
ejecta and to a SN remnant for Newtonian ejecta. Or, bright transient emission, dubbed a 
'mini SN' (Li & Paczyhski 1998; Rosswog & Ramirez- Ruiz 2002), is produced by radioactive 
elements that are synthesized during the rapid decompression of very dense and neutron rich 
material that is ejected during a NS-NS or a NS-BH merger (see, e.g., Rosswog et al. 1999). 
Our upper limit at t ~ 22.8 days constrains mainly the former mechanism, in particular 
the amount of energy in ejecta with r > 1.3(E 51 /noY^ 8 , and suggests that the total energy 
in a relativistic outflow is significantly smaller in GRB 050509B than that in typical long 



3 A higher £?k,iso is possible for a very low external density. For n <~ 10 6 cm 3 , i?k,iso ~ 10 51 erg for 
z = 0.225 and - 10 52 - 53 erg for z w 3 (Bloom et al. 2005; Lee, Ramirez-Ruiz & Granot 2005). This would, 
however, require -E 7j i S o <C -Ek.iso, i.e., a very inefficient prompt emission (compared to -E 7) i so ^ £-k,iso for long 
GRBs), and would not naturally reproduce the fact that Lx/E y ^ so for GRB 050509B is similar to that for 
long GRBs. 
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GRBs. 4 Our upper limit at t — 1.85 days provides more stringent constraints on the latter 
mechanism, in which the emission is expected to peak around the optical-UV range within 
a day or so (up to a few days) with a semi-thermal spectrum (Li & Paczyhski 1998). The 
'mini SN' emission is mainly concentrated in a very narrow energy range (i.e., the optical) 
during (and near) the peak; therefore, the X-ray emission could have been easily missed by 
the Chandra observation of GRB 050509B at t ~ 2.5 days. 

Using the simplified model of Li & Paczyhski (1998), the optical flux from a 'mini-SN' as- 
sociated with GRB 050509B should have been a factor of ~ 10 3 (//10- 3 )(M/0.01M Q ) 1 / 2 (3t;/c) 1 / 2 
higher than our upper limit at t = 1.85 days, where M and v are the mass and velocity 
of the ejected material, and / is the fraction of its rest energy that goes into radioactive 
decay of the elements with a decay timescale comparable to the time it takes the expand- 
ing ejecta to become optically thin, t T . For a kinetic energy of 10 51 _E 5 i erg, where E 51 = 
(M/0MM & )(3v/c) 2 ~ 1, varying M and v within a reasonable range (0.003 < M/M Q < 1 
and 0.03 < v/c < 0.5) would not change the optical luminosity by more than one order of 
magnitude. A larger uncertainty is the value of /, which reflects the amount of radioactive 
material synthesized in the accompanying NS-NS wind. From the above simple arguments 
we derive an approximate upper limit of / < 10~ 5 . We note here that the most efficient 
conversion of nuclear energy to the observable luminosity is provided by the elements with a 
decay timescale comparable to t T . In reality, there is likely to be a large number of nuclides 
with a very broad range of decay timescales. Our limits thus place interesting constraints 
on the abundances and the lifetimes of the radioactive nuclides that form in the rapid de- 
compression of nuclear-density matter - they should be either very short or very long when 
compared to t T so that radioactivity is inefficient in generating a large luminosity. 

The above arguments suggest that either the intrinsic energy in the outflow from 
GRB 050509B was <C 10 51 erg, or alternatively, and arguably more likely, that it was close 
to the canonical value of ~ 10 51 erg but most of this energy was in sub-relativistic ejecta 
with a very small radioactive component during the optically thick expansion phase. 5 The 
latter is very different from long/soft GRBs which typically have ~ 10 51 erg in highly rela- 
tivistic ejecta with r > 100. More short GRB afterglows are needed to establish whether 
the energetics of GRB 050509B is representative of the bulk of the short/hard GRB class. 

Finally, our observations may place constraints on other possible models for short GRBs 



4 One possible caveat is the dependence of the afterglow brightness on the density of the burst environment 
(see Bloom et al. 2005); since the possible counterpart location on Gl spans a large range of densities, we 
have not folded this dependence in our conclusions. 

5 Such sub-relativistic velocities could be the result of a significant entrainment of baryons into the outflow. 
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(e.g., Ramirez- Ruiz 2004; Tan, Matzner & McKee 2001; Rosswog, Ramirez- Ruiz & Davies 
2003; Shibata & Sekiguchi 2003; Blackman & Yi 1998; Rosswog et al. 2004; Davies, Levan & 
King 2005). The strict optical upper limits derived in this Letter, argue that these scenarios 
are only feasible if the transport of the energy is in the form of subrelativistic ejecta with 
little or almost no radioactivity, or in any other form of delayed energy input such as provided 
by a pulsar or by later mass ejection by a central source. 

We thank Elena Pian and Alberto Castro-Tirado for comments and the ESO Paranal 
Science Operations staff, in particular Chris Lidman and Andreas O. Jaunsen for efficiently 
conducting the observations reported here. The Dark Cosmology Centre is supported by the 
DNRF. The authors acknowledge benefits from collaboration within the EC FP5 Research 
Training Network "Gamma- Ray Bursts - An Enigma and a Tool" . 
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Fig. 1. — Left: First epoch V image (3.92 days after burst) showing the putative host 
elliptical galaxy Gl and several faint galaxies in the XRT error circle. Middle: Same, after 
subtraction of a fit to Gl. The cross marks the location of the center of Gl. North of it is 
the new detected source which may be a foreground or background source or a companion 
to Gl. Right: Difference between last (22.86 days after burst) and first epoch V images. 
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20 




Time since GRB (days) 

Fig. 2. — The upper limits (arrows) on variable sources (in the difference images) inside the 
GRB 050509B XRT error circle at the epochs given in Table 1 compared to the lightcurves 
of different SNe redshifted to z — 0.225. Solid curves indicate Type Ic SNe, dashed curves 
Type la SNe. Thick solid curve: The hypernova SN 1998bw accompanying the long GRB 
980425. Thin solid curve: The faint Ic supernova SN 19941. Thick dashed curve: A typical 
Type la SN (stretch = 1). Thin dashed curve: A faint Type la SN similar to SN 199 lbg. A 
Galactic extinction of E(B - V) = 0.019 mag (Schlegel et al. 1998) towards GRB 050509B 
has been assumed. 
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Table 1. Log of observations 



Date 


Phase 


Band 


Exp. time 


Seeing 


Lim. mag 


(UT) 


(Days past GRB) 




(s) 


(arcsec) 


(mag) 


050511.02 


1.85 


R 


2700 


0.9 


26.6 


050513.09 


3.92 


V 


2700 


0.9 


27.5 


050517.10 


7.93 


V 


1800 


0.7 


25.0 


050522.05 


11.88 


V 


1800 


1.0 


24.2 


050601.00 


22.83 


R 


2700 


0.9 


26.7 


050601.03 


22.86 


V 


2700 


1.0 


27.5 



Note. — The quoted 3a limiting magnitudes are measured in the field in a 
2xFWHM aperture. The limiting magnitudes become progressively smaller 
towards the center of Gl. 



